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Abstract — In this paper we investigate the optimal relay 
location in the sense of maximizing suitably defined coverage 
region for MIMO relay channel. We consider the general 
Rayleigh fading case and assume that the channel state 
information is only available at the receivers (CSIR), which is 
an important practical case in applications such as cooperative 
vehicular communications. In order to overcome the mathe- 
matical difficulty regarding determination of the optimal relay 
location, we provide two analytical solutions, and show that it is 
possible to determine the optimal relay location (for a desired 
transmission rate) at which the coverage region is maximum. 
Monte Carlo simulations confirm the validity of the analytical 
results. Numerical results indicate that using multiple antennas 
increases coverage region for a fixed transmission rate, and also 
increases the transmission rate linearly for a fixed coverage. 

Index Terms — Optimal relay location; coverage region; MIMO 
relay channel; desired transmission rate. 

I. Introduction 

The relay channel is the most basic structural unit 
in wireless networks and relaying strategy can increase 
channel throughput, coverage region, facilitate information 
transmission and realizes some of the gains of multiple- 
antenna systems by single-antenna terminals. The relay 
channel, since its introduction by van der Muelen [1], has been 
extensively studied [2] -[6]. In their seminal work, Cover and El 
Gamal [2] presented a capacity upper bound and achievability 
strategies for the relay channel. Although the channel capacity 
is still unknown, the authors in [6] unified most of known 
capacity theorems into one capacity theorem, which potentially 
may be applicable to a more general class of relay channels. 

Using multiple antennas can be considered as an effective 
technique to combat fading which also can result in a 
significant increase in channel throughput. Nonetheless, the 
capacity gain obtained from this technique heavily depends 
on the amount of instantaneous CSI available at the receivers 
and transmitters [7]. Considering different scenarios of CSI, 
the capacity for point-to-point MIMO channel in [8] and for 
MIMO relay channel in [9]-[ll] have been analyzed. 

The existing results for MIMO relay channels mainly 
focus on maximizing capacity bounds for fixed locations 
of the channel nodes. However, in most wireless networks, 
the optimal location of the relay and coverage region is of 
practical interest due to the mobility of the destination (e.g. 




Fig. 1. MIMO relay channel. 



a mobile station in a cellular network). The authors in [12] 
have investigated single-antenna Gaussian relay channel with 
the new objective of maximizing coverage for a desired 
transmission rate, and by considering decode-and-forward 
(DF) and compress-and-forward (CF) strategies for the relay 
channel, optimized the relay location as a design parameter. In 
this paper, we define the concept of coverage and investigate 
the coverage region in a more general and practical case, i.e., 
MIMO fading relay channel with only CSIR. More precisely, 
our goal is to determine the optimal relay location (for a 
desired transmission rate) at which the coverage region is 
maximum. Thus, considering the coverage definition and the 
necessary condition for applying DF strategy, we express the 
desired transmission rate in terms of the optimal relay location 
and derive two expressions with the help of which we can 
determine the optimal relay location. 

The rest of the paper is organized as follows. We start with 
the channel model in section II. In section HI, by defining the 
channel path-loss coefficients, we review the capacity bounds 
of MIMO relay channel in an appropriate form for the case 
when the channel entries are i.i.d. Rayleigh fading and only 
CSIR is available. The concept of coverage and evaluation of 
optimal relay location are provided in section IV. Numerical 
results are presented in section V, and section VI contains our 
conclusion. 

Throughout this paper, scalars are represented by lowercase 
letters, vectors are denoted by lowercase boldface letters and 
uppercase boldface letters are used for matrices. Superscripts 
T and H denote the transpose and conjugate transpose, respec- 
tively. The operator £[.) stands for the expectation and log(.) 
denotes base-2 logarithm. I m is the M x M identity matrix, 
and we use tr(A) and det(A) to denote the trace and the 
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Fig. 2. Network geometry. 



determinant of the matrix A, respectively. The distribution of 
a circularly symmetric complex Gaussian vector x with mean 
fi and covariance matrix Q is denoted by x ~ CAf(fi, Q), 
and 2 ~ C7V(0, 1) is a circularly symmetric complex Gaussian 
random variable, where its real and imaginary parts are zero 
mean i.i.d. Gaussian random variables, each with variance 1/2, 
i.e., JV(0,l/2). 

II. Channel Model 

Consider the MIMO relay channel in Fig. 1. The source 
transmits a message X\ over the channel. Let us consider y 2 
as the received signal at the relay; based on prior received 
signals, the relay then transmits a message X2 that is intended 
to facilitate the transmission between the source and the 
destination. We assume the source has Mi transmit antennas, 
and the destination has N 3 receive antennas. Also, suppose 
that the relay is equipped with M 2 and N 2 antennas for 
transmitting and receiving, respectively. Thus, the received 
signals at the destination and the relay can be expressed as 



V 3 = H sl xi + H 32 x 2 + n 3 



y 2 = H 2 ix 1 



n 2 



(1) 



(2) 



where 



• x\, x-2 are Mi x 1 and M 2 X 1 transmitted signals of the 
source and the relay, respectively; and have zero mean, 
i.e., £[xi] = 0, S[x 2 ] = 0; 

• 2/3 > 2/2 are JV3 x 1 and N 2 x 1 received signals at the 
destination and the relay, respectively; 

• H31, if 21, H 32 are N 3 x Mi, N 2 x Mi, and N 3 x M 2 
channel gain matrices as shown in Fig. 1, and modeled 
as independent (fiat) fading processes; 

• n 3 , n-2, are respective independent ^3 x 1 and N 2 x 1 
zero-mean circularly symmetric complex Gaussian noise 
vectors at the destination and the relay, with distribution 
C7V(0, 1 n 3 ) and CAf(0, In 2 ), which are independent of 
the transmit signals. 

Considering M = Mi + M 2 , the M x M joint transmit 
covariance matrix of the zero-mean source and relay transmit 
signals can be described by 



Q 



Qn Q12 
Q21 Q22 



(3) 



where = £[xiX^],i, j = 1,2 is the covariance matrix 
between the input signals Xi and Xj. Note that Q is a 
Hermitian matrix [8]. Throughout this paper, we suppose that 
the relay works in full-duplex mode (the relay can transmit 
and receive in the same frequency band at the same time), 
because it can be considered as a performance upper bound 
for half-duplex systems [11]. 

Now consider the network geometry depicted in Fig. 2. 
In this configuration the source and the relay are located at 
di = (0, 0) and d 2 = (u 2 , 0), respectively, and the destination 
is located at d 3 = (u 3 ,v 3 ). Also, we suppose that 77 is the 
distance-based path-loss power attenuation exponent. Thus, the 
channel matrices can be represented as follows 



H 



31 



a-H w3 i, if 21 — b.H w2 i, H 32 — c.H w32 (4) 



where the channel path-loss coefficients defined as: 



a = (u§ + vD^W, b 4 



A — v/ 2 A 



C A ( (ua _ U3 )2 + v 2)-u/4 > 



the entries of H w3 \, H w2 \ and H w32 are i.i.d. £/V(0, 1), and 
each block use of the channel corresponds to an independent 
realization of channel matrices. 

III. Capacity Bounds of MIMO Relay Channel 

In this section, we review the capacity upper bound and 
lower bound for MIMO relay channel, and present the results 
obtained for the case of only CSIR. 

A. Cut-Set ( CS) Upper Bound 

Using the "max-flow min-cut" theorem, Cover and El 
Gamal [2] showed that the capacity of the full-duplex relay 
channel in terms of the channel mutual information is upper 
bounded by 

C C s = max min{I{xi]y 2 ,y 3 \x 2 ), I{x u x 2l y 3 )} (5) 

where the maximization is with respect to the joint 
distribution of the source and relay signals. Considering 
Xi ~ CAf(0, Qa), i = 1,2, where Q ti is the covariance matrix 
of Xi, the mutual information expressions in (5) can be 
expressed for the MIMO relay channel as [9] 



C C s = max 



Ci = logdet(J 



N 



C 2 = logdet(J 



min(Ci,C 2 ) (6) 

£=1,2 

H 1 Q 1{2 H?) (7) 

+ H 2 QH%) (8) 



where Hi 



H, 



H 



31 



H 



32 



N = N 2 + N 3 



H 3 \ 
H 2 \ 

and Q 1 | 2 = £[x\X± \x 2 ] = Q n — Q\ 2 Q 22 Q 2 \ is the 
conditional covariance matrix and given by Schur complement 
of Q 22 in Q [13]. The optimal distribution p(x\,x 2 ) in (5) is 
Gaussian [2], and consequently the maximization of (6) would 
be with respect to three covariance matrices Q n , Q 22 , and 

Ql2- 

When the channel matrices are random (due to fading) and 
the CSI is only known at the receivers, the optimal joint 



transmit covariance matrix Q in (6) is diagonal. Using Jensen's 
inequality, the authors in [14] showed that the equal power 
allocation is the optimal solution, i.e., 



Q 



22 



P2 T 



Q 



12 







(9) 



where Q 12 — refers to the independence between the source 
and the relay signals. Thus, the CS upper bound for the MIMO 
relay channel with only CSIR can be expressed as 



C§ s =mm(C?,C?) 



Cf = £ 



log det(J 



N 



Mi 



(10) 



(11) 



logdet(/jv 3 + H 2 



Ml' 1 Mi 



P2 

M 2 



H H 2 ) 



(12) 

where the superscript R stands for "Rayleigh", and the 
expectations are taken over the channel gains matrices. 

B. DF Achievable Rate 

The capacity of the full-duplex relay channel is lower 
bounded by the DF achievable rate [2] given as 

Rdf = max min {/(an; y 2 \x 2 ), I(xi, x 2 \ y 3 )} (13) 

p(a:i,a:2) 

where the optimal distribution is again Gaussian and the 
maximization should be done over the joint distribution of the 
source and relay signals. In this strategy the relay first decodes 
the received signal from the source, and then re-encodes it 
before forwarding it to the destination. Considering (4) and 
using the same approach for evaluation of the CS upper bound, 
it is easy to show that the capacity of MIMO relay channel 
with only CSIR is lower bounded by 



R% F = unn(C*,C*) 



£ 



log det(J 



\ j -t- -j^j--b 2 -H w21 H H 



w2\> 



(14) 



(15) 



where C 2 is the same as in (12). Note that (15) is the ergodic 
capacity of MIMO channel between the source and the relay, 
and we will use it in the following section to find the optimal 
relay location in the sense of maximizing coverage region. 

IV. Main Results 

In this section, considering a desired transmission rate, 
we define the concept of coverage in MIMO relay channel. 
Next, we investigate the coverage region and the optimal relay 
location based on: 1) an exact expression for the ergodic 
capacity of MIMO channels, and 2) an approximation in the 
high-SNR regime. Finally, we examine the impact of using 
multiple antennas on the optimal relay location. 



A. Coverage Definition 

Our definition of coverage region has a close relation to the 
concept of outage capacity [15], and we consider it for MIMO 
relay channel as a geographic region at which a rate of at least 
R > is guaranteed, i.e., 



A(u 2 ) = {d 3 :C(u 2 ,d 3 )>R} 



(16) 



where R denotes the desired transmission rate in bps/Hz, 
C(u 2 , d 3 ) is the channel capacity when there is a fixed distance 
u 2 between the source and the relay, and the destination is 
located at d 3 = (u 3 ,v 3 ). 

Since the capacity of MIMO relay channel is still an open 
problem in general, it can be inferred from the condition 
in (16) that the lower bound of the channel (DF achievable 
rate) should be larger than R. On the other hand, for applying 
DF strategy it is necessary for the relay to have the ability of 
decoding the information which is transmitted to it from the 
source. Thus, in order for the R to be achievable by the source- 
relay channel, considering ( 1 5), the desired transmission rate 
R can be defined as follows 



R^£ 



log det(J 



N 2 



■jj-ti w21 ti w21/ 



(17) 



where p = Pi.d*~ v is the effective SNR, and d* is a boundary 
distance which determines whether the DF strategy can be 
applied or not, since the condition u 2 < d* guarantees that 
the relay is still able to decode the transmitted signal from 
the source. We define d* as the optimal relay location (for 
a desired transmission rate R) at which the coverage region 
is maximum. Our intuition for this finding is based on the 
fact that when the relay locates at distances smaller than d* 
{u 2 < d*) the coverage region decreases, while for distances 
larger than d* (u 2 > d*) the DF strategy can not be used 
anymore in accordance with our definition of coverage in (16). 
Therefore, u 2 — d* is the best choice for relay location in the 
sense of maximizing coverage. 

B. Desired Transmission Rate Analysis 

Note that (17) involves the expectation operator that 
generally admits no explicit solution. In order to obtain a 
theoretical expression between the desired transmission rate R 
and the optimal relay location d*, in what follows, we evaluate 
R by using two analytical approaches. 

1) Exact Expression for Desired Transmission Rate: As 
stated before, the R in (17) is the capacity of source-relay 
MIMO channel, thus it can be expressed by the exact ergodic 
capacity of uncorrelated Rayleigh fading MIMO channel with 
only CSIR as [16, Theorem 9] 



R 



1 



ln2.r m (m).r m (n) 



5^det(*(i)) (18) 



1=1 



where m — min(Mi, N 2 ), n = max (Mi, N 2 ), and the 
modified multivariate Gamma function T t (k) and the auxiliary 
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Fig. 3. The relation between R and d* when M\ = JV2 = 2. 




Fig. 4. The relation between d* and d^ when R = R±, and Mi = = 4. 



to x to matrix ^P(Z) are defined as follows 

t 

r t (k)^l[T(k-p + i) 



*(0 = 




Mi ' ' 



i = Z 



(19) 



(20) 



where 1 < i,j < m, s = n — m + i + j — 1, T(fc + 1) = kl, 
and the auxiliary function JF{e, h,g) is defined by 

9 i 

•F(e,fc,$) ± h'S.e^ .J2(h.er.T(-g+ Pl —) (21) 

p=i 

with incomplete Gamma function T(—k,l/x) for positive 
integers k defined as in [16, eq. (90)]. 

2) 1: Now we restrict our attention to the high-SNR regime, 
and derive the following theorem. 

Theorem 1: In the high-SNR regime, the relation between 
the desired transmission rate R and the optimal relay location 
d* can be approximated by the following expression 



R w to. log 



p. exp 



Mi 



m n— p 

lyy- 

ln 2 ^ ^ q 

p=i 9 =i y 



(22) 



where —"0(1) ~ 0.577215 is the Euler-Mascheroni constant. 

Proof: For large p we can ignore the identity matrix /at 2 
in (17), and simplify it as follows 



R 



£ 



£ 



logdet(^H w21 H^ 21 ) 



logdet(^-H^ 21 H w21 ) 



Mi > N 2 



M x < N 2 . 



Considering the case Mi > N 2 , R can be rewritten as 



R - N ^Wi + i-2 £ 



\ndct(H w2 iH" 21 ) 



(23) 



(24) 



Note that H w2 iH w2 i is an i.i.d. Wishart matrix with M 2 
degree of freedom and covariance matrix I^ 2 . Thus, the 
expectation in (24) can be expressed as (see [17, Theorem 
2.11]) 

N2-1 



£ 



lndet(H w2 iH% 2 i)\ = ^ i>{M 2 - p) (25) 

p=0 



where ip(r) is Euler's digamma function, which for natural r 
can be expressed as [18] 



q=l 



Thus, (25) can be rewritten as 

£ 



(26) 



N 2 M 2 -p 



\ndet(H w21 H% 21 )\ = N 2 .^) +J2 E ~ ^ 

p=l q=l q 

and applying (27) to (24), the result follows. Note that by using 
the determinant identity det(7 + AB) = dct (I + BA) and 
following similar steps in (24)-(27), the desired transmission 
rate can be evaluated for the case Mi < N 2 . 

C. The Impact of Using Multiple Antennas on the Coverage 

In this section, considering the high-SNR regime, we 
investigate the effect of using multiple antennas on the 
coverage region of MIMO relay channel in Rayleigh fading 
environment. Thus, we express d* in terms of the optimal relay 
location in single-antenna relay channel. 

Using (22), the desired transmission rate for the single- 
antenna relay channel (Mi = N 2 = 1), can be written as 



Ri wlogtPi.dp'-exp^W) 



(28) 



where d\ denotes the optimal relay location in this channel. 
Considering (22) and (28) for a fixed transmission rate 
(R = Ri), d* can be expressed as follows 



d* = 



(Pi.exp^ 1 >) m - 1 .df .2 J 



MY' 



where A is defined as 



In 2 ^ ^ q 

p=l q =l 



(29) 



(30) 



Remark: Note that d* in (29) has an inverse relationship with 
the number of transmit antennas, and thus we can conclude 
that receiver diversity provides wider coverage region than 
transmitter diversity. 
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Fig. 5. Coverage region when 112 = 0.95, and N = 2. 
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Fig. 6. Coverage region when 112 = 1.05, and N = 2. 

V. Numerical Results 

In this section we present the numerical results for the 
coverage region of MIMO relay channel. In our simulations 
we assume that 7] = 3.52, Mi = N 2 = M 2 = N 3 = N, and 
P\= P 2 = 10 dB. Fig. 3 depicts the desired transmission rate 
R for different values of d*. The region below each curves 
represents the region where DF strategy can be applied, while 
the region above contains the points in which this strategy 
can not be used anymore. It can be seen from the figure that 
there is a perfect match between Monte Carlo simulation and 
the analytical expression in (18), and also the approximation 
in (22) is more accurate at high-SNR than at low-SNR. 

In Fig. 4 the effect of using multiple antennas on the optimal 
relay location is investigated. This figure shows the ratio d* Jd\ 
for different values of transmit power versus the optimal relay 
location for single-antenna relay channel. It can be inferred 
from the figure that as the distance between the source and 
the relay gets larger, the effect of using multiple antennas 
decreases. 

For desired transmission rate R = 5.5, the optimal relay 
location is d* = 1 [see (18)]. Fig. 5 shows the coverage region 
for U2 = 0.95 while in Fig. 6 it is assumed that u 2 = 1.05. 
For performance comparison, these figures also include the 
coverage region of the "no-relay" cases. Fig. 5 indicates that 
if u 2 < d* the capacity bounds converge but the coverage 
decreases. On the other hand, if u 2 > d*, it can be inferred 
form the Fig. 6 that although the CS upper bound may extend, 
but in this case the DF strategy can not be used anymore. 
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Fig. 7. Coverage region when R = 4, and 112 = 0.775. 




No. of antennas 



Fig. 8. Transmission rate vs. number of antennas for a fixed coverage. 

Fig. 7 shows the growth of coverage region for a fixed 
transmission rate with increasing the number of antennas, and 
Fig. 8 depicts the growth of transmission rate for a fixed 
coverage. It can be observed that for a fixed coverage region 
there is a linear relationship between the transmission rate and 
the number of antennas. 

VI. Conclusion 

In this paper we investigated the optimal relay location in 
the sense of maximizing coverage region for MIMO relay 
channel. Considering the case when the channel matrix entries 
are i.i.d. Rayleigh fading and the CSI is only available at the 
receivers, we defined the coverage region for MIMO relay 
channel, and derived an exact analytic expression with the help 
of which we could determine the optimal relay location (for a 
desired transmission rate) at which the coverage region is max- 
imum. Also, an approximation is presented for the high-SNR 
regime. Numerical results confirm the accuracy of our analysis, 
and show that using multiple antennas increases coverage 
region for a fixed transmission rate, and also increases the 
transmission rate almost linearly for a fixed coverage region. 
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